Astronomy <fe Astrophysics manuscript no. bbhyp 


February 2, 2008 


(DOI: will be inserted by hand later) 





Hyperon softening of the EOS of dense matter and the spin 
evolution of isolated neutron stars 

J.L. Zdunik 1 , P. Haensel 1,2 , E. Gourgoulhon 2 , and M. Bejger 1 



CO 

o 
o 

CN 

> 
O 

o 

CN 



> 
O 



cn 
o 

On; 

6 ■ 
a: 

C3 



N. Copernicus Astronomical Center, Polish Academy of Sciences, Bartycka 18, PL-00-716 Warszawa, Poland 
2 LUTH, UMR 8102 du CNRS, Observatoire de Paris, F-92195 Meudon Cedex, France 

jlz(5camk.edu.pl, haensel@camk.edu.pl, Eric.Gourgoulhon@obspm.fr, bejger@camk.edu.pl 

Received / Accepted 

Abstract. The effect of the hyperon softening of the equation of state (EOS) of dense matter on the spin evolution 
of isolated neutron stars is studied for a broad set of hyperonic EOSs. We use a multidomain 2-D code based 
on a spectral method, and show how important the precision of solving the equations of stationary motion is for 
the stability analysis. For some EOSs, the hyperon softening leads to spin-up by the angular momentum loss, in 
the form of the back bending phenomenon, for a rather broad range of stellar baryon mass. We show that large 
segments of the evolutionary tracks exhibiting the back bending behaviour in the moment-of-inertia - rotation- 
frequency plane, are actually unstable and therefore not astrophysically relevant. We show also that during the 
spin-up - angular-momentum-loss epoch, an isolated neutron star (e.g., a radio pulsar) can lose a sizable part of its 
initial angular momentum without significantly changing its rotation period. We propose also simple arguments 
and criteria allowing to connect the presence of a back bending epoch with the mass-radius relations and the 
stiffness and/or softness of the nucleon and hyperon EOSs of the neutron star core. 
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1. Introduction 

Some theories of dense matter predict a softening of the 
equation of state (EOS) at densities exceeding normal 
nuclear density po — 2.7 x 10 14 g cm~ 3 , implied by a 
phase transition to a new "exotic" phase. Several "ex- 
otic" high-density phases were proposed in the past, e.g., 
deconfined quark plasma, and pion-condensed or kaon- 
condensed hadronic matter. The softening of the EOS 
could be due to a transition into a pure "exotic" phase, 
or to a mixture of an "exotic" phase with normal phase of 
dense matter. It has been suggested that the softening of 
the EOS due to a phase transition could lead to character- 
istic back bending phenomenon in the timing of spinning- 
down pulsars (Glendenning et al. 1997) or produce char- 
acteristic period clustering in spinning- up neutron stars in 
low-mass X-ray binaries (Glendenning 2001, Glendenning 
and Weber 2002). An interesting conclusion of these pa- 
pers was that both back bending in spinning down pulsar 
and spin clustering in accreting millisecond pulsars stars is 
an evidence of a phase transition taking place at the cen- 
ter of a spinning-down pulsar or a spinning-up accreting 
neutron star. 

As we show in the present paper, using several EOS of 
dense matter and very precise code for the calculation of 



the rotating stellar configurations, observation of the back 
bending in the timing behavior of isolated pulsars, or of 
the period clustering, is not an unambiguous evidence for 
an "exotic" phase in dense matter. Namely, these phe- 
nomena can also be implied by the presence of hyperons 
in dense matter, a feature which is in no way "exotic" and 
which was predicted more than forty years ago (Cameron 
1959, Salpeter 1960, Ambartsumyan and Saakyan 1960). 
A possibility of "spin-up by the angular momentum loss" 
for a normal sequence (baryon mass smaller than the max- 
imum baryon mass of static configurations) of spinning- 
down neutron stars with hyperonic cores was previously 
noted by Balberg et al. (1999). However, as we will demon- 
strate in the present paper, a complete study of back bend- 
ing requires a very precise "exact" code for calculating 
stationary configurations of rotating stars and a simulta- 
neous careful checking of secular stability of these config- 
urations. Such conditions were typically not satisfied in 
previous calculations. 

In previous works the back bending phenomenon has 
been considered as a feature of the I(fi) dependence, 
where / is the moment of inertia of the star and O is the 
angular frequency of rotation ( Glendenning et al.(1997)[ 



Cheng et al.(2002) and Spyrou & Stergioulas (2002) \. In 
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the present paper we clarify some statements about the 
back bending for rotating neutron stars, and we formu- 
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late some rules which are useful for searching for the back 
bending in rotational stellar sequences. As we demonstrate 
using high- precision evolutionary sequences, the whole 
analysis should be performed using different pair of vari- 
ables: total stellar angular momentum J versus ft instead 
of 1(0,). High precision is particularly important because 
it is needed to reliably check the secular stability of ro- 
tating configurations; only the stable ones are interesting 
and observationally relevant. 

The paper is organized in the following way. Softening 
of the EOS due to the presence of hyperons is discussed 
in Sect. |21 The method allowing for a high precision of the 
2-D calculations of the equilibrium configurations of rotat- 
ing neutron stars is briefly described in Sect. [31 Different 
formulations of the stability criteria for rotating configu- 
rations are briefly summarized in Sect. In Sect. |3]we 
propose a method of checking for the occurrence of the 
back bending phenomenon by inspecting the baryon-mass 
- equatorial-radius relations at fixed values of rotation fre- 
quencies. We apply this method to several EOSs with a 
hyperon softening. The interplay between the back bend- 
ing and stability is discussed in Sect. [fJl where we study 
neutron-star evolution tracks in the angular-momentum - 
rotation-frequency plane. In Sect. [S| we study the depen- 
dence of the back bending phenomenon on the EOS of the 
hyperonic matter. Final Sect. contains discussion of our 
results, including their possible observational aspects, and 
ends with concluding remarks. 
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Fig. 1. Equations of state of hyperonic matter calculated 
by Balberg & Gal (1997) and used in this paper. Our 
notation is a close analogue of that introduced by Balberg 
et al. (1999). Our labels Nl, N1H1, and N1H2 are their 
EoSl N, EoSl NAS, and EoSl NAES, respectively. Our 
N2H1 and N2H2 are their EoS2 NAS and EoS2 NASS, 
respectively. For further explanations see the text. 



2. Equation of state with hyperons 

Possible presence of hyperons in dense neutron-star matter 
is mainly the consequence of the Pauli principle for neu- 
trons and electrons, which at a sufficiently high baryon 
density can make a replacement of high-energy neutrons 
and electrons by more massive but slower hyperons en- 
ergetically favorable. Hyperon species H is then present 
above a certain threshold pn, which is determined by a 
condition involving the change of energy of dense matter, 
due to addition of a single hyperon, at a fixed pressure, 
H° H . The threshold density for H is the lowest density at 
which the equality 



(1) 



is satisfied. Here, p n and p e are the chemical potentials 
of neutrons and electrons (which include rest energies), 
and qn is the hyperon charge in units of e. The threshold 
condition given by Eq. 0) shows that because of the high 
value of fi e , hyperons with negative charge are strongly 
preferred over the positively charged ones. It also explains, 
why S _ , and not the lowest mass hyperon, A , is usually 
the first hyperon to appear in dense matter. It should be 
stressed, however, that the hyperon-nucleon interactions 
which contribute significantly to p? H , are poorly known, 
and this implies uncertainty not only in the values of pu 
but also in the order of the hyperon appearance. Let us 
notice that for some models of dense matter hyperons do 
not appear in neutron stars at all, because the lowest ph 



is larger than the maximum density in neutron stars (see, 
e.g., Pandharipande & Garde 1972). 

The hyperon softening is particularly well visualized 
by the behavior of the adiabatic index 



r 



n h dP 



(2) 



At each threshold density p = pn the function T(p) suf- 
fers a drop. For some models, typical values of T ~ 2 — 3 
characteristic of the nucleonic matter at p ~ (2 — 4)po 
can drop even down to T ~ 1 (see, e.g., Balberg et al. 
1999). For other models the drop is not so dramatic but 
still sizable (see, e.g., Haensel et al. 2002). The softening 
of the EOS by hyperonization of matter can be also visu- 
alized by comparing the values of the maximum allowable 
mass for non-rotating neutron stars, M max , for the EOS 
without hyperons, referred hereafter as the N EOS, and 
those involving nucleons and hyperons (NH). Typically, 
allowing for the presence of hyperons lowers the value of 
M max by (0.3-0.6) M Q (see, e.g., Haensel 2003). The pres- 
ence of hyperons leads to a very characteristic flattening of 
the mass-radius and mass-central density plot for neutron 
star, with a knee taking place just after the threshold for 
the first hyperon (Glendenning 1985, Balberg et al. 1999). 
This feature will be important in the context of the back 
bending phenomenon in rotating neutron stars. 

In present paper we use the EOSs calculated by 
Balberg and Gal (1997). These EOSs are presented in Fig. 
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n They are based on phcnomcnological effective baryonic 
matter energy functionals, resulting from effective baryon- 
baryon interactions whose parameters are adjusted to re- 
produce the empirical properties of nuclear matter as well 
as the basic experimental features of the hyperonic inter- 
actions. Two models of nucleon matter, based on an ef- 
fective nucleon-nucleon (NN) interaction, constructed by 
Balberg and Gal (1997) lead to the EOSs Nl and N2 of 
dense matter without hyperons (they correspond to their 
models EoSl N and EoS2 N, in the notation of Balberg et 
al. 1999). The model Nl corresponds to the incompress- 
ibility of nuclear matter at saturation Ko = 240 MeV (i.e., 
is close to the experimental value of this parameter) . The 
model N2 gives K = 320 MeV, and leads to a signifi- 
cantly stiffer EOS than the Nl one. Balberg & Gal con- 
structed also two models representing effective nucleon- 
hyperon (NH) and hyperon-hyperon (HH) interactions in 
baryonic matter. We will denote them by HI and H2; they 
correspond to the AS and ASS effective hyperon interac- 
tion models of Balberg & Gal (1997) in the notation of 
Balberg et al. (1999). Model H2 leads to a stronger soft- 
ening of the EOS due to the presence of hyperons than the 
HI one. The EOS of baryonic matter are obtained by com- 
bining NN, NH, anf HH effective interactions. In this way 
one obtains four EOSs of baryonic matter with hyperons. 
We will denote them by N1H1, N1H2, N2H1, and N2H2; 
they correspond to EoSl NAS, EoSl NAES, EoS2 NAS 
and EoS2 NAES in the notation of Balberg et al. (1999), 
respectively. For the sake of completeness, we considered 
also five EOSs of Glendenning (1985); see Sects. El and 03 

3. Numerical calculations 

The neutron star models have been computed in 
full general relativity by solving the Einstein equa- 
tions for stationary axi-symmetric spacetime (see e.g. 
Bonazzola et al. (1993) or Gourgoulhon et al. (1999) for 



approach is automatic fulfillment of the first law 
of thermodynamics (the Gibbs-Duhcm relation) (see 
Nozawa et al. (1998)|). 



the complete set of partial differential equations to be 
integrated). The numerical computations have been per- 
formed via the Lorene/Codes/Rot_star/rotstar code 
from LORENE (http://www.lorene.obspm.fr/). This 
C++ code implements a multi-domain spectral method 



introduced in |Bonazzola et al. (1998) A description of the 
code can be found Gourgoulhon et al. (1999) For the pur- 
pose of the present work, we have employed two domains 
to describe the neutron star interior, making use of the 
adaptive coordinates to set the boundary between the in- 
nermost domain and the outer one at the transition sur- 
face to hyperon matter. In this way, the density field is 
smooth in each domain and the spectral method results 
in a high accuracy. This accuracy has been checked by 
evaluating the GRV2 and GRV3 virial error indicators (see 
Nozawa et al. (1998) I, which showed a relative error lower 
than ~ lO" 15 . 

The physical parameters resulting from the equation 
of state (pressure, energy density, number density) are 
obtained by the Hermite interpolation ||Swesty (1996) 
Nozawa et al. (1998)). The important feature of this 



4. Back bending phenomenon, I and J vs. Q, and 
stability of rotating configurations 

The term back bending comes from nuclear physics 
(see, e.g., Ring & Shuck 1980). Nuclei can be excited by a 
projectile to a state of a rapid rotation corresponding to 
a large angular momentum quantum number J and ex- 
citation energy E(J). The nuclear angular momentum is 
measured in the units of fi. The eigenvalues of the operator 
of the square of the angular momentum are J(J+1). For 
J ^> 1 one can approximate this eigenvalue by a "clas- 
sical value" J 2 . Within the quasiclassical approximation 
(in which J can be treated as a continuous quantity) one 
can phenomenologically define an "angular frequency" by 
to = dE/dJ. The nuclear moment of inertia T is then 
found by fitting the rotational E(J) spectra. In the quasi- 
classical approach, X is a function of to, which can be rep- 
resented by a curve in the X -co plane: along this curve, J 
increases. In the standard case there is a one to one cor- 
respondence between X and to, and both X(co) and lu(X) 
are increasing functions of their arguments. However, for 
some nuclei (e.g., 158 Er see Ring & Schuck 1980) to reaches 
a maximum at some value of J = J\ and then decreases 
(back bends) to reach minimum at some larger value of 
J = Ji- At J = Ji and J = J 2 the derivatives duj/dj 
and dui/dX vanish, and for J\ < J < J 2 the curve to(X) 
is "back bending" which corresponds to to which is de- 
creasing with increasing I (see Fig. 3.4 in Ring & Schuck 
1980). 

Neutron star can be treated as a huge atomic nucleus. 
However, such a star is a macroscopic classical object con- 
taining some ~ 10 57 baryons, to be compared with at 
most ~ 200 — 250 nucleons in rapidly rotating nuclei. In 
nuclei, all high-angular-momentum states are the excited 
ones and therefore unstable. In the case of rotating neu- 
tron stars unstable states are not interesting, their lifetime 
being too short to observe them. Therefore in the case 
of neutron star we have to check whether a given state 
of stationary rotation is stable, because instability would 
make it astrophysically irrelevant. In this context, we find 
it convenient to discuss the back bending phenomenon by 
studying the dependence of the total angular momentum 
of the star J versus frequency of rotation / = oj/2tt. We 
differ in this choice from the previous work, in which the 
dependence /(/) was studied. The reason of this choice is 
the following: J is a well defined quantity describing in- 
stantaneous state of rotating relativistic star and the evo- 
lution of rotating star can be easily calculated under some 
assumptions about the change of J. In what follows, we re- 
strict ourselves to evolutionary tracks of isolated neutron 
stars for which baryon mass Mb = const. 

The moment of inertia / is usually defined as J /VI (see, 
e.g., Stergioulas 2003). This is the definition used in the 
previous papers on the back bending phenomenon in ro- 



4 



J.L. Zdunik et al.: Hyperons and back bending in neutron stars 




f [IIz] 

Fig. 2. The difference between the back bending curves 
J(f) and /(/)• The spin-evolution track is calculated for 
M B = 2.15 M Q and N2H1 EOS. Upper panel presents 
the dependence J(/). The thicker segment AB represents 
stable back bending evolution, while dotted segment to 
the right of the minimum at B consists of configurations 
unstable with respect to axisymmetric perturbation. The 
lower panel represents the /(/) spin evolution track, con- 
sidered in the previous work; the back bending behavior 
occurs apparently along the whole AC branch. Actually, 
the BC segment of the /(/) is astrophysically irrelevant, 
because configurations to the right of B are unstable. 



fating neutron stars. However, such / does not describe the 
response of the star to the change of J or Q, and therefore is 
not useful for checking the stability of rotating configura- 
tion. In order to obtain, e.g., the spin down of a star due 
to the angular momentum decrease dJ one should have 
defined / = dJ/dQ. Only in the slow rotation limit, where 
only terms linear in are conserved, both definitions of 
the moment of inertia coincide. 

Total stellar angular momentum J is not only a quan- 
tity with a strict physical meaning in general relativity. It 
also allows us to study the stability of rotating configu- 
rations with respect to axially symmetric perturbations. 
The point of the change in stability within a family of 
rotating configurations (from stable to unstable or vice 
versa) . corresponds to the extremum of M or Mb at fixed 
J IjFriedman et al. 1988|> : 



dx I 







J —const 



dAh 
dx 



0, (3) 



J— const 



where x is the first of the two parameters which 
parametrize (label) stationary rotating stellar configura- 
tions (the second one being J), for example x = p c or 



Fig. 3. Total baryon mass Mb vs circumferential equato- 
rial radius R cq for stationary configurations rotating at a 
fixed frequency, for the N2H1 EOS. The curves are labeled 
by the rotational frequency / = Q/2n [Hz]. The dotted 
line corresponds to the onset of instability with respect 
to axi-symmetric (quasi-radial) perturbations. All curves 
terminate on the large-radius side at the mass-shedding 
(Keplerian) configurations. The maximum baryon mass 
of non-rotating stars, Mg*^ ax , is marked by a dash-dotted 
line. 



x = P c . Equivalently, the onset of instability can be de- 
termined by the condition 



dJ 
dx 







Mr —const 



dJ_ 

dx 







(4) 



M— const 



The crucial difference between the dependencies J(f) 
and /(/) in the back-bending region is presented in Fig. 
it will be discussed in more detail in Sect. El The ap- 
parently "back bending branch" in lower panel, as defined 
in the previous work (Glendenning et al. 1997), consists 
mostly of configurations unstable with respect to axial 
perturbations. The BC segment is astrophysically irrel- 
evant, and the B point correspond to an unstable termi- 
nation of the back bending branch (see Sect.0). 

5. Softening of the EOS and the 

baryon-mass-equatorial-radius relation at fixed 
rotation frequency 

Rotation modifies the relation between the baryon mass 
(Mb) and the circumferential equatorial radius (i? C q) for 
equilibrium stellar configurations. The baryon mass Mb 
plays a special role, because it remains constant during 
the evolution of solitary pulsars. In the present section we 
point out specific features of the MB(i? cq ) curves, which 
signal the presence of the back bending (BB) phenomenon 
in rotating neutron stars. Our calculations were performed 
for all the equations of state with hyperons presented 
in Balberg et al.(1999) However, we present detailed re- 
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R eq [km] 

Fig. 4. The enlarged region of the back-bending phe- 
nomenon (corresponding to the box in Fig. OJ) . For fixed 
rotational frequency /i,i n fx — 880 Hz and /2,infx — 1320 Hz 
the Mb(R cci ) dependence has an inflexion point (cor- 

( 1 infx) 

responding to the masses Mg' ' = 1.91 M and 
Afg 2 ' mfx ' ) = 2.17 Mq, respectively) resulting in a region 
where the curve is nearly flat. For / S [/i,i n fx, /2,infx] there 
exists a local minimum of Mb- Dashed curves correspond 
to a fixed total angular momentum J. The configurations 
close to the local maxima are obviously stable (for a fixed 
J, Mb is monotonic). The evolution of an isolated star 
which is losing its angular momentum is represented by 
the motion along a horizontal line from right to the left 
(decreasing J). The loss of J in the back bending regime 
is associated with a spin up of the star. 




i i i i i i i i i i i i i i i i i i i i i i_ 

8 10 12 14 16 18 

R eq [km] 

Fig. 5. Same as in Fig. [3 but for the N1H1 EOS. 



suits (figures) only for those EOSs for which the BB phe- 
nomenon is strongly pronounced. 

In Fig.[3]we show the Me(i? C q) curves for stars rotating 
at a fixed rotation frequency / = Sl/2%, calculated for the 
N2H1 EOS. An enlarged view of a particularly interesting 




R eq [km] 

Fig. 6. The enlarged region of the back-bending phe- 
nomenon (corresponding to the box in Fig. EJ. 



rectangular region of the Mb — Req plane is shown in Fig. 

5.1. Signature of BB: minimum of the baryon mass at 
fixed frequency 

As we will show, the BB phenomenon is strictly con- 
nected with the existence of a local minimum of Mb in the 
Mb (-Req) plot at a fixed /. The softening of the EOS due 
the hyperonization leads to the flattening of the Mb (.Req) 
and Mb(/0c) curves (for the case of non-rotating stars, see 
Balberg et al. 1999). This effect of flattening grows with 
increasing rotational frequency and for / > /l.infx rota- 
tion may even produce a local minimum of MB(R C q) and 
M(R cq ). The curve M B (R cq ) [or M B (p c )} at / = / Unfx 
has a very specific property. Namely, for this curve first 
and second derivatives of Mb with respect to the central 
density p c vanish at some p c — /9 c ,mfx, i-e., the curve has 
there a point of inflexion which corresponds to baryon 
mass Mb m . Our numerical calculations for the N2H1 
and N1H1 EOS give / Unfx ~ 880 Hz and f XMx ~ 860 Hz 
respectively (see Figs. 0] and EJ- For / > /i,infx the curve 
MB(-Req), Mb(Pc)) ■ ■ ■ j exhibit a local minimum. We find 
that the presence of this local minimum is an indication 

( 1 infx) 

that for Mb > M B ' the rotational evolution of neu- 
tron star with J < exhibits a BB phenomenon in the 
vicinity of / ~ /i,i n f x - This is clearly seen in Figs. 0] and 
[HI where we show an enlargement of the vicinity of the 
inflection point, together with J = const, lines. As we see, 
in this case there exist a range of Mb where the decrease 
of J leads to the increase of the angular frequency which 
is exactly equivalent to back bending. The fragment of 
the curve for which M decreases as a function of p c does 
not necessarily correspond to the instability region - the 
decrease of Mb at a fixed / does not imply the decrease 
at a fixed J. It is only the latter condition which indi- 
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cates the instability with respect to small axi-symmetric 
perturbations. 1 

In Figs. 0] and El we draw also three horizontal lines 
corresponding to fixed values of the total baryon num- 
ber. A rotating star losing its angular momentum moves 
along horizontal line from the right to the left. The bottom 
lines correspond to the mass Mb = Mj' lnfx ' ~ 1.91 M© 
(~ 1.81 M ) for the N2H1 (N1H1) EOS, at which the 
curves for / = A.infx — 880 Hz (860 Hz) have a point 
of inflexion. The top horizontal line corresponds to the 
different situation in which the curve Mft(R eq ) has an in- 
flexion point at a higher frequency, namely at / = 1320 
Hz and 1100 Hz, for the N2H1 and N1H1 EOSs, respec- 
tively. For baryon masses larger than the mass at this 
inflexion point, the angular momentum loss does not lead 
to the decrease of angular frequency before the onset of 
instability is reached, the star is all the time accelerat- 
ing. The value of this limiting masses are Mb — 2.17 M© 
and 1.90 M Q , respectively. The intermediate horizontal 
line corresponds to the maximum mass of the non- rotating 
stars M B = M B *^ ax = 2.05 M Q and 1.824 M Q . 

5.2. Acceleration or slowing down close to the 
Keplerian limit 

The interesting difference between the cases N2H1 and 
N1H1 concerns the behavior of the rotating star as it starts 
losing J at the Keplerian frequency /k- The question is 
whether the star is then slowing down or spinning up. The 
actual behavior can be easily deduced from the shapes of 
the M B (i? eq ) [or Mb(/O c )] curves for a fixed frequency close 
to the Keplerian limit. If Mb is increasing as we move in 
the Mb — i? q plane away from the Keplerian configuration 
[(dMs/d/3 c )K > 0], the star is slowing down as it loses the 
angular momentum. This effect results in the S-shape of 
the J(f) dependence in the back-bending case. Otherwise 
[(dAfB/dp c )K < 0] the isolated star is spinning up when 
evolving from the Keplerian configuration with J < 0. The 
limiting case corresponds to the condition (dMB/dp c )K = 
[or (dMB/di? cq )K = 0] and is represented in Fig. 0]by 
the curve with / = 970 Hz (the baryon mass at this point 
M |K,flat) = 2 Q1 For the Nim EOg guch an effect 

appears for supramassive stars (/ ~ 1030 Hz, M B K ' flat ' ) = 
1.96 M©, see Fig. El 

6. Back bending and stability: analysis in the 
J — f plane 

The search for the BB phenomenon with simultaneous 
testing of the stability of rotating configurations can be 
most conveniently carried out by plotting, at a fixed Mb , 
the stellar angular momentum J versus rotation frequency 
/. Let us start with the N2H1 EOS where the BB is the 

1 It may be shown that the configurations which realize the 
extrema of the Mb and M at fixed J coincide, see Harrison 
et al. 1965 for the static case, and Friedman et al. (1988) for 
uniformly rotating configurations. 



most pronounced. Several curves J = «/(/) at selected val- 
ues of Mb, calculated for this EOS, are shown in Fig. 
These curves represent the evolution of an isolated pulsar 
of a given baryon mass Mb, as it loses its angular mo- 
mentum due to radiation of electromagnetic waves. Along 
each curve, the central density p c increases monotonically 
when one moves downward. For stable configurations J 
is a monotonic function along this path. Any minimum 
indicates the onset of the instability with respect to axi- 
symmetric perturbations. Putting it differently, for stable 
configurations each value of J corresponds to one and only 
one value of /. The BB manifests itself as a stable segment 
of the J(f) curve with dJ/d/ < 0. 

It should be mentioned that there is an important dif- 
ference between the information one can get from the 
analysis of the J(f) and the usually used /(/) curves. 
Although the /(/) can have segments with dl/df < cor- 
responding to a back bending in a nuclear physics sense 
(Glendenning et al. 1997), the /(/) dependence cannot tell 
us whether a seemingly "back bending branch" contains 
configurations which are stable. One can have a minimum 
in J(f) on a "back bending" segment of /(/) (see Fig. 
EJ) ; such a possibility was already mentioned by Spyrou & 
Stergioulas (2002). 

The final fate of the rotating star as its angular mo- 
mentum decreases depends on M B . If M B > M B ^* ax , the 
star is supramassive and eventually collapses into a black 
hole. 

For the N1H1 EOS the BB phenomenon is less pro- 
nounced than for the N2H1 one. The spin-evolution tracks 
at Mb = const are shown in Figs. E| [</(/) curves] andlTfH 
[zoomed J(f) curves in the BB region for normal con- 
figurations]. The zoomed fragment in Fig. 1101 shows how 
narrow, compared with the N2H1 case, is the range of 
baryon masses for normal configurations where BB with 
stable termination occurs. 

Using the method developed in Sect.[21 one can readily 
formulate the following criterion. If M B 1 ' mfx ^ > M B *^ ax 
then the back bending is possible only for supramassive 
configurations. In the opposite case, the back bending can 
occur also for rotating stars with Mb < M B ^ ax . 

7. Two types of back bending 

A look at Fig. shows that there are actually two types 
of BB for rotating neutron stars. For isolated neutron star 
J < and therefore on a stable J(f) branch J 2 < J\ 
implies t 2 = t(J 2 ) > t\ — t(Jx). Assume then that a 
BB epoch starts at t\ and terminates at t 2 . For the first 
type of BB, occurring for the baryon mass 2.25 M Q and 
2.2 M© for the N2H1 EOS (Fig. Q, and for 2 M© and 
1.95 M© for the N1H1 EOS Fig.©, the point at t 2 corre- 
sponds to an unstable rotating configuration. Such a BB 
with an unstable termination will be denoted as BB ut . 
This type of BB occurs above some value of M B equal to 
M B 2 ' lni *\ An isolated neutron star enters the BB regime 
but it will not reappear in the normal spin-down regime 
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Fig. 7. Angular momentum of the star J versus rota- 
tion frequency /, for the N2H1 EOS. Each curve corre- 
sponds to a fixed Mb, whose value in M Q is displayed. 
Along each curve, the central density increases down- 
wards. The dotted segments correspond to configurations 
which are unstable with respect to axi-symmetric pertur- 
bations whereas the thick lines correspond to spin-up by 
angular momentum loss. An enlarged view of the rect- 
angular region within which the back bending occurs is 
shown in Fig. [5] 



with dJ/d/ > 0. However, the situation changes as de- 
creasing Mb crosses the value M^' mtx \ Then the BB seg- 
ment splits into two segments. The segment with higher J 
(corresponding to an earlier epoch) terminates at a stable 
configuration, and will be denoted by BB s t. Then follows 
an epoch of normal spin down with d J/df > 0, followed by 
a BB ut terminated by an instability. With decreasing Mb, 
the BB ut epoch becomes shorter and shorter, and becomes 
infinitesimally small as one reaches the maximum mass 
of nonrotating configuration Afg" ffl . For baryon masses 

(1 infxl 

larger than M B ' we are dealing with a "BB episode" 
in an otherwise normal neutron star rotational evolution. 

The case of BB ut can be quite easily discussed on the 
basis of the Me(i?eq) dependence (for example Fig. |3J. 
The BB ut phenomenon corresponds to a specific location 
of the instability line (dotted line in Fig. 0} with respect 
to the maxima of the functions MB(i? C q) at fixed / (if 
they exist). Because these maxima are to the right of the 
instability line (i.e., in the stable region) BB ut always ap- 
pears. In the other words, a star with fixed Mb, approach- 
ing the instability point, has to eventually spin up. This 
conclusion follows immediately from a Lemma formulated 
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Fig. 8. Enlargement of Fig. 




500 1000 1500 



f [llz] 

Fig. 9. Same as FigEJbut for N1H1 EOS. The dash-dotted 
line is the J(f) curve for the Nl EOS (i.e., not allowing 
for the presence of the hyperons). 

by Friedman et al. (1988). 2 These authors considered a 
two-parameter family of uniformly rotating stars with a 
one-parameter EOS P = P(p). In general, a continuous 
sequence of rotating configurations can be labeled by a 

2 We note a misprint in this Lemma text as printed in 
Friedman et al. (1988) , who use dots to denote a derivative with 
respect to A: the dot over the right bracket in (fl J + fiN) ^ 
is missing. 
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Fig. 10. Angular momentum of the star as a function of 
frequency for the N1H1 EOS, in the region where the back- 
bending phenomenon with stable termination (BB st , see 
the text) occurs. Each curve corresponds to a fixed value 
of Mb (in M©). Along each curve, the central density in- 
creases downwards. Segments along which the angular mo- 
mentum loss is associated with a spin up are indicated by 
a thick solid line. All configuration on the S-shaped thick 
segments are stable. The dash-dotted line describes a J(f) 
trajectory for the EOS of matter in which hyperons are 
not allowed (Nl EOS, see the text). This curve illustrates 
the dramatic effect of the hyperon softening of the EOS 
on the pulsar spin-evolution. 

parameter A, so that along this sequence all stellar pa- 
rameters are functions of A. In order to avoid confusion 
with time derivative of a stellar quantity Q, we will de- 
note a derivative of Q with respect to A by Q 1 = dQ/dA. 
According to Friedman et al. Lemma, the unstable region 
corresponds to the part of the sequence </(/) for which 
J'f > (since along these sequences M B = 0), or equiv- 
alently, dJ/d/ > 0, which means that just before the in- 
stability is reached one has dJ/d/ < 0, i.e., the spin up 
by the angular momentum loss. 

Concluding, BB ut is a feature of any equation of state. 
However the significance of this effect depends on the stiff- 
ness of the matter. In contrast BB st exists only when Mb 
has a local minimum for fixed rotational frequency. 

The impact of the hyperon EOS softening on the char- 
acter and significance of the BB phenomenon can be 
clearly appreciated by comparing Figs. with Fig. ITU 
obtained for the (nucleonic) Nl EOS without hyperons. 
For the nucleonic EOS, BB is present only for supramas- 
sive configuration and it is significant only for relatively 
high masses (although strictly speaking it is present for 
all masses larger than - for masses Mb < 2.6 M© 
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Fig. 11. Angular momentum of a star with fixed baryon 
mass Mb as a function of frequency for the Nl EOS. 
Segments constraining unstable configurations are dotted, 
those corresponding to the spin-up by the angular momen- 
tum loss are indicated by a thick solid line. 
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Fig. 12. Model N2H2. Baryon-mass versus circumferen- 
tial equatorial radius for fixed frequency. The minimum for 
fixed Mb appears approximately at the same frequency at 
which the decreasing part for large R eq (close to Keplerian 
configurations) disappears, i.e. this limiting case corre- 
spond to the nearly flat curve (Mb = const) up to the 
Keplerian configuration. There is no maximum for fixed 
frequency. 



the BB ut region is so narrow that it is not visible in Fig. 
ITn . see Cook et al. (1994) and Salgado et al. (1994). In 
contrast, in the case of the EOS with hyperons, the BB 
is present for configurations close to the Keplerian one for 
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Fig. 13. Model N2H2. The enlarged region marked by the 
rectangular box in Fig. El For / ~ 920 Hz there exist re- 
gion where the curve is nearly flat (Mb does not change) 
up to the Keplerian configurations. For higher / there ex- 
ist a minimum of Mb . The evolution of the isolated star 
which is losing its angular momentum is represented by 
the motion along horizontal line from right to the left (de- 
creasing J). For M-q > 1-82 M angular momentum loss 
leads to spin up. For explanation of the "way features" of 
most curves see the text. 

very broad range of masses. The BB epoch for the nucle- 
onic EOS terminates always by an instability. 

8. The EOS-dependence of the back bending 
phenomenon 

The specific form of the BB phenomenon depends on the 
degree of softening of the EOS above the hyperon thresh- 
old. We considered a rather large set of the hyperonic 
matter, including four models of Balberg & Gal (1997) 
described in Sect. and the five models of Glendenning 
(1985). Rapid rotation and spin-evolution for the N1H1 
and N2H1 EOSs of Balberg & Gal (1997) and of an EOS 
of Glendenning (1985), were previously studied, using dif- 
ferent numerical method of solving of the 2-D equations 
of hydrostatic equilibrium, by Balberg et al. (1999). 

We obtain the most pronounced BB behavior for the 
N2H1 EOS, which was not studiedby Balberg et al. (1999). 
In this model the nucleonic EOS is very stiff, and the 
hyperon softening is not so strong as for the EOS with 
the H2 hyperonic sector. Such features are optimal for 
getting a strong BB effect for a large interval of stel- 
lar baryon mass. In particular, we get BB st for a rela- 
tively large interval of normal rotating configurations with 
M B = (1.91 - 2.05) M Q , which extends up to 2.17 M 
if supramassive rotating models are also included. At 
Mb — 2 Mq the star which starts its spin evolution at 
/ ~ /k accelerates its rotation by ~ 50 Hz during the 
first epoch of its evolution when nearly half of its initial 
angular momentum is lost ! Implications of such a sub- 



stantial "back bending with stable termination" for the 
pulsar timing interpretation will be discussed in Sect. [UJ 

The N1H1 EOS has significantly softer nucleon com- 
ponent, and more important flattening of the MB(i? C q) 
curve. As for the N2H1 EOS, the BB st phenomenon oc- 
curs not only for the supramassive configurations but 
also for the normal ones, with M B < Mg ^ ax (Fig. EJ). 
However, the range of the baryon masses of the normal 

configuration for which the BB st occurs is very narrow, 
M a,i„fx) = Lgl Mq < Mb < M stat^ = L823 Mq For 

such baryon masses BB st takes place at / > /i jm fx — 
850 Hz. The BB st itself for this EOS is due to the high 
degree of flatness of the hyperon section of the Mef-Rcq) 
plot. However, due to this flatness the values of Afg 1,lnfx ^ 
and M^^ax are also very close to each other. All in all, 
even including the supramassive configurations, the range 
of baryon masses for which the BB st occurs is not large, 
1.81 Mq < M B < 1.9 M Q . 

The situation changes if we consider the N2H2 EOS 
(see Sect.EJ), where the hyperon softening is stronger. Let 
us repeat the analysis applied previously to N1H1 EOS 
in Sect. 03 The MB(i? cq ) curves at several values of / are 
shown in Fig.fR] In contrast to the case of the N1H1 EOS, 
pictured in Fig. there is no curve with a visible second 
(large-radius) maximum. A zoom of the relevant baryon- 
mass - equatorial-radius rectangle is displayed in Fig. 1131 
A marginal second maximum appears at / ~ 920 Hz. 
However, this second maximum is not a robust one. It 
can be detected only if the precision of the calculation is 
sufficiently high. This is to be contrasted with the N1H1 
case, where the second maxima on the Mb — R eq curves 
are well pronounced, see Fig. |SJ Let us notice that N2H2 
is the EOS for which the BB phenomenon for normal con- 
figurations was detected by Balberg et al. (1999). 

The main features of the dependence of the back 
bending phenomenon on the EOS can be summarized on 
the basis of the Table containing relevant parameters 
for considered EOSs. Let us start with four characteris- 
tic baryon masses M^' infx) , M^' infx) , M B K ' flat) , M B ^ ax . 
Within this set, the ordering according to the baryon- 
mass value turns out to depend only on the pure nu- 
cleon EOS. The sequence for N1H1 and N1H2 EOSs is: 

M (l,infx) < ^ < M (2,infx) < ^K.flat) for 

the N2H1 and N2H2 EOSs we have M^ 1,infx) < A4 K ' flat) < 
max < Mg ,mfx ' ) . As can be seen the main difference is 
a rather low value of M^' Rlit) for the N1H1 and N1H2 
models. In the case of the N2 model of the nucleon com- 
ponent, rotation changes the properties of the star close 
to the Keplerian frequency much more effectively than for 
the Nl one. 

The nucleon N2 EOS is stiffer than the Nl one. As 
a result, the stars with nucleon envelopes based on the 
N2 model are more extended than the Nl ones, hav- 
ing larger radius at the same mass. On the other hand, 
the density profiles in the nucleon envelopes of the N2H 
(H=H1,H2) stars have smaller radial gradients, and there- 
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Table 1. Main parameters relevant for the spin evolution 
of isolated neutron stars with hyperons. The labels "infx" 
mark the configurations for which the curve AfB(J? eq )/ 
has a point of inflexion. For /i^nfx < / < /2,infx (corre- 
sponding masses M B Unfx) and M^ 2 ' infx) ) BB st exists. For 
Mb > M§ m&x we have also BB ut , which is the only BB 
phenomenon for Mb > M^' lni *\ For / > /ic.fiat and 
Mb > Afg K,flat \ a star which is initially in the mass- 
shedding (Keplerian) state, spins-up as it looses angular 
momentum. The characteristic S-shape of J(/) is then im- 
possible. 



Parameter 


N1H1 


N2H1 


N1H2 


N2H2 


m^sumq] 


1.83 


2.05 


1.73 


1.95 


/l.infx [Hz] 


860 


880 


880 


910 


A4 1,fa&) [M Q ] 


1.808 


1.91 


1.73 


1.81 


/K.fiat [Hz] 


1030 


970 


1050 


930 


M (,K,flat) j M@ j 


1.96 


2.01 


1.87 


1.82 


/2,infx [Hz] 


1100 


1320 


1080 


1380 


M,< 2 ' infx) [M ] 


1.9 


2.17 


1.79 


2.06 



fore they play more important role during rotation. For 
example, the N2 envelopes have significantly larger mass 
than the Nl ones. Consequently, a smaller rotational fre- 
quency (/K.flat = 930, 970 Hz for H2, HI softening, respec- 
tively), than in the case of the Nl model, is needed to make 
the effects of rotation so important that the maximum of 
M close to Keplerian frequency disappears. This effect 
can be easily seen in Table where not only the N1H fre- 
quencies /K,flat are significantly larger than the N2H ones, 

but also the mass M B K ' flat ^ is significantly larger. Not only 
it exceeds the maximum mass for the nonrotating stars, 
M" a , but it is even larger than M^' mix \ In this case 

the M^' in determines the mass above which only spin 
up by angular momentum loss is possible. 

As the Nl EOS is softer than the N2 one, the addi- 
tional softening by hyperons leads almost immediately to 
the maximum mass of nonrotating configurations; the hy- 
peron segment beyond the "hyperon knee" is very flat. 
The back bending phenomenon BB st is present between 
the two inflexion points, and in general is not connected 
with the maximum mass of nonrotating stars. 

The MB(Rcq) curves displayed in Fig. IT31 deserve an 
additional comment, referring to the precision of the 2- 
D calculations. Nearly all curves for / = 850 — 1100 Hz 
(except for two curves for / ~ 920 Hz) were obtained 
in numerical calculations in which the innermost zone 
boundary is not adjusted to the surface of the hyperon 
threshold. Consequently, these curves exhibit "wavelets" 
which result from an insufficient precision of the numerical 
calculations. In contrast, the two curves for / ~ 920 Hz 
are calculated with the innermost zone boundary at the 
hyperon-softening threshold, which enables a much higher 
precision. 



The division of the stellar interior into two zones is 
particularly effective in the case of a strong change in the 
EOS P(p) at some density. In our case we have a rather 
stiff EOS below the threshold density for the hyperon ap- 
pearance, and a soft EOS for dense matter with hyperons. 
Therefore, we put the inner zone boundary at the thresh- 
old of hyperon appearance, where the adiabatic index of 
the EOS suffers a significant drop. This allows us for very 
accurate calculation by the spectral method also in the re- 
gion close to the hyperon threshold. Although in our case 
we do not encounter a density jump at the zone boundary, 
our method can be also used for an EOS with a density 
jump due to a first-order phase transition. 

9. Discussion and conclusion 

The presence of hyperons in neutron-star cores can 
strongly affect the spin evolution of a solitary neutron 
star (e.g., an isolated pulsar). As we have shown, such 
a neutron star can live a long epoch of "spin down by 
the angular momentum loss" , and this could occur for a 
broad range of baryon mass of neutron star. The epoch 
with P < despite J < could terminate by an instabil- 
ity or by a stable continuous transition to a "standard" 
spin-down epoch. 

We paid particular attention to the (P < 0, J < 
0, Mb = const) epoch with a stable termination. It is 
represented by an S-shaped segment of the spin-evolution 
track in the J — / and I — f planes, and was baptized 
"back bending phenomenon" in the previous literature. 
Various regimes of the spin evolution were shown to be 
correlated with the behaviour of the MB(-R e q)/ curves at 
fixed /. In particular, we pointed out importance of the 
location of the inflection points of the MB(R C q) curves for 
the existence of the back bending phenomenon. We were 
also inclined to leave the name of back bending only to the 
evolution-track segments with a stable termination (BB st 
in our terminology). 

Epochs with back bending for normal rotating config- 
urations were found for two of the four EOS of Balberg 
& Gal (1997). On the other hand, we found that the back 
bending phenomenon for normal sequences was absent for 
five hyperonic EOSs of Glendenning (1985). This illus- 
trates the uncertainties in the hyperonic EOS, stemming 
for a high degree of ignorance concerning the nucleon- 
hyperon and particularly hyperon-hyperon interactions in 
dense matter. 

Throughout this paper we stressed that the back bend- 
ing in the I — f plane, considered up to now in the liter- 
ature, should be accompanied by stability analysis; only 
stable back-bending configurations are astrophysically rel- 
evant. We performed such a stability analysis, and we 
found that very often dominant back-bending segments of 
the /(/)m b are unstable with respect to the axisymmetric 
perturbations, and therefore do not exist in the Universe. 

As it has been mentioned by Spyrou & Stergioulas 
(2002), in such kind of calculations it is extremely im- 
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portant to assure the thermodynamical consistency of the 
EOS (the first law of thermodynamics has to be strictly 
fulfilled), ft is well known that a rough treatment of this 
condition can lead to an inaccurate determination of the 
maximum mass and stability conditions (for example the 
configurations corresponding to the maximum of M and 
Mb do not coincide - this is an evidence for a lack of the 
thermodynamical consistency in the EOS). 

The precision of the code is also very important and 
in our case the proper division of the star into two com- 
putational domains (at the threshold for appearance of 
hyperons) allows us to obtain high precision results with- 
out a large increase of the number of grid points. 

An isolated pulsar, born in a SN II explosion, could 
have an interesting and nonstandard past due to a hy- 
peron softening of the EOS. As we showed, such a pul- 
sar could lose some half of its initial angular momen- 
tum without changing much its rotation period. Therefore, 
if one observes a rapid pulsar with a characteristic age 
t psr significantly longer than the age of a supernova 
remnant where this pulsar is born, tsnr, this might be 
due to some back bending episode, resulting from a (hy- 
peron? phase-transition?) softening of the EOS of its 
core. Such a possibility of explaining a seeming contra- 
diction between tsnr and tpsr has been already no- 
ticed, in the context of a mixed-phase EOS softening, by 
Spyrou & Stergioulas (2002) Clearly, the observational 
pulsar-timing signatures of the EOS softening due to the 
hyperon or phase-transition softening of the EOS deserve 
further studies, and we are planning to continue such stud- 
ies using our high-precision 2-D code. 

Another interesting consequence of the hyperon soften- 
ing of the EOS could be a "period clustering" of rotating 
neutron stars powered by accretion in the long living low- 
mass X-ray binaries. This problem is now being studied. 
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